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Deep radio polarization studies, together with detailed 
MHD modeling, can independently verify those conclu- 
sions, that are based on Hi observations and dynamical 
models. We performed deep polarized radio-continuum ob- 
servations of the Virgo cluster spiral galaxy NGC 4654 
with the Effelsberg 100m telescope at 8.35 GHz and the 
VLA at 4.85 GHz. Detailed 3D MHD simulations were 
made to determine the large-scale magnetic field and the 
emission distribution of the polarized radio continuum in 
the model, during the galaxy evolution within the clus- 
ter environment. This direct comparison between the ob- 
served and simulated polarized radio continuum emission 
corroborates the earlier results, that the galaxy had a re- 
cent rapid close encounter with NGC 4639 and is undergo- 
ing weak ram pressure by the intracluster medium. This 



^v^j ■ combination of deep radio polarization studies and de- 



tailed MHD modeling thus gives us unique insight into 
the interactions of a galaxy with its cluster environment. 
It represents a diagnostic tool that is complementary to 
deep Hi observations. 
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Abstract. A recent comparison between deep VLA Hi observations and dynamical models of the Virgo cluster 
spiral galaxy NGC 4654 has shown that only a model involving a combination of a tidal interaction and ram 
pressure can reproduce the data. 



1. Introduction 

The bright Virgo cluster spiral galaxy NGC 4654 is located 
at the periphery of the cluster at a projected distance of 
3.4° = 1 Mpc 1 from the cluster center (M87). Its two- 
armed spiral structure is asymmetric with a shorter spiral 
arm to the northwest and a larger wider spiral arm to 
the southeast. Its Hi distribution shows a sharp edge to 
the northwest and a large, very low column density tail 
to the southeast (Phookun & Mundv ITMSjl . While the 
northwestern part of the Hi velocity field corresponds to a 
rising rotation curve, rotation is almost absent in the very 
low column-density tail. The Ha image (Koopmann et al. 
2001J shows enhanced star formation at the northwestern 
edge of the galactic disk. 

Another spiral galaxy, NGC 4639, is located only 18' = 
90 kpc in projection from NGC 4654 and might be consid- 
ered to be a companion galaxy of NGC 4654. NGC 4639 
has a ring-like Hi distribution. Since it is located at the 
edge of the field observed by Phookun & Mundy l|1995[) 
and Cayatte et al. I|1994f) . it is not clear if its gas distribu- 
tion and kinematics are perturbed. While Warmels ( 1988) 
observed extended gas in the east of NGC 4639, Cayatte 
et al. (1994J) found atomic gas to the southwest. 

Two kinds of interactions can in principle be respon- 
sible for the observed asymmetries and perturbations of 
NGC 4654: (i) a gravitational interaction with a mas- 
sive galaxy or (ii) ram pressure due to the galaxy's rapid 
motion within the hot and tenuous intracluster medium. 
In order to investigate the recent evolutionary history of 
NGC 4654, Vollmer (2003) made numerical simulations 



1 We use a distance of 17 Mpc to the Virgo cluster 



including a gravitational interaction, ram pressure, and a 
combination of both, concluding that 

— a past ram-pressure stripping event where the galaxy 
has already passed the cluster center 800 Myr ago and 
is now emerging from the cluster core gives rise to a 
low surface-density tail in the southeast of the galaxy 
center. This tail shows clear signs of rotation; 

— a gravitational interaction between NGC 4654 and 
NGC 4639 cannot produce a low surface density tail; 

— only a gravitational together with a relatively small 
amount of ram pressure leads to the observed low 
surface-density tail where rotation is almost absent. 

The results of Vollmer PIUS}) suggest that NGC 4654 
had a tidal interaction with another massive galaxy in the 
recent past (~500 Myr ago). A model where the position 
and radial velocity of the perturbing massive galaxy are 
reasonably close to those of NGC 4639 can explain the per- 
turbed stellar distribution and the gas distribution of the 
inner disk. Therefore, Vollmer pTOlj) identified NGC 4639 
as the perturbing galaxy. 

A caveat of the gravitational interaction scenario is 
that NIR Tully-Fisher measurements yield a line-of-sight 
distance between NGC 4654 and NGC 4639 of more than 
5 Mpc (Gavazzi et al. lTMty . While NGC 4654 is located 
1-2 Mpc in front of the cluster core, NGC 4639 is lo- 
cated far behind (~ 5± 1 Mpc). Observations of Cepheids 
give independent distance estimates for a number of Virgo 
cluster galaxies (Gibson et al. 2000, Freedman et al. 2001). 
Five galaxies are in both, the Gavazzi et al. ( 1999) and the 
Freedman et al. H2001|l samples. For 4 galaxies (NGC 4321, 
NGC 4535, NGC 4536, and NGC 4536), the difference be- 
tween the distance moduli obtained from the two methods 
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are small (|A/Xo| < 0.2). However, this difference is much 
larger for NGC 4639 (|Ayu | = 0.4). The Cepheid distance 
estimates place NGC 4639 closer, but still a few Mpc be- 
hind the center of the Virgo cluster (M87). Thus these 
measurements corroborate the large NIR Tully Fisher dis- 
tance of NGC 4639. 

Since the preferred scenario of Vollmer et al. (2003) 
includes ram pressure of p lam — 200 cm _3 (kms _1 )~ 1 , a 
galaxy velocity of ~ 1000 kms -1 implies an intraclus- 
ter medium (ICM) density of nicM ~ 10~ 4 cm~ 3 . At 
the distance of NGC 4639, the ICM density according 
to Schindler et al. ljH?9"9"|) is lower than 10~ 5 cm" 3 . This 
suggests that either (i) NGC 4654 had a gravitational 
interaction with another massive object, (ii) the Tully 
Fisher method underestimates the line-of-sight position 
of NGC 4654, (i.e. NGC 4654 is too bright with respect 
to its rotation velocity), or (hi) the ICM distribution is 
elongated along the line-of-sight, similar to the galaxy dis- 
tribution (Gavazzi et al. 11999(1 . In the case of scenario (i) 
the perturber has to be a massive dark galaxy (see, e.g., 
Minchin et al. I2005[) since there is no other visible com- 
panion than NGC 4639. 

In this article we test the viability of the preferred 
scenario of Vollmer (2003J). We used polarized radio- 
continuum observations of NGC 4654, together with an 
MHD model, to put further constraints on the models. 
The radio continuum emission is due to relativistic elec- 
trons that gyrate around the galactic magnetic field (syn- 
chrotron emission). If the polarized radio-continuum emis- 
sion (called PRCE hereafter) is free of Faraday rota- 
tion (typically at frequencies above ~ 4 GHz), it traces 
the ordered large-scale magnetic field directly. This mag- 
netic field is extremely sensitive to shear and compres- 
sion motions, which are not easily detected in 2D velocity 
fields. Using polarized radio-continuum emission at 6 cm, 
Vollmer et al. $2004) could unambiguously identify the 
ram-pressure compressed region in the Virgo cluster spi- 
ral galaxy NGC 4522. We were able to reproduce the Hi 
distribution, velocity field, and the distribution of polar- 
ized radio-continuum emission with our model simulations 
(Vollmer et al. 2006). 

In order to understand the behavior of the large- 
scale magnetic field and thus the polarized radio- 
continuum emission during a ram-pressure stripping 
event, Otmianowska-Mazur & Vollmer (2003J solved the 
induction equation on the velocity fields provided by a dy- 
namical model. They showed that a prominent maximum 
of polarized radio emission is observed during compression 
and in a state where the galaxy has passed the cluster 
center, the ram-pressure decreases, and the gas pushed by 
ram pressure falls back onto the galactic disk. We applied 
the method described in Otmianowska-Mazur & Vollmer 
(2003J) to NGC 4654; i.e., we solve the induction equation 
on the velocity fields of the simulations of Vollmer ( 2003 ) 
and compare the results to our polarized radio-continuum 
observations. 

The Effelsberg and VLA radio continuum observations 
are presented in Sect.[3 In Sect.Elwe explain the 3D MHD 



model and how we calculated the model of polarized radio- 
continuum emission maps. The results of our model calcu- 
lations are shown in Sect. 0] followed by direct comparison 
between our observations and simulations (Sect.[5J). These 
results are discussed in Sect. and we give our conclusions 
in Sect. □ 

2. Radio observations of NGC 4654 

Radio polarimetric observations of NGC 4654 were per- 
formed at 8.35 GHz using the single-horn receiver in the 
secondary focus of the 100-m Effelsberg radio telescope 2 . 
The galaxy has been scanned alternatively in R.A. and 
Dec. Four data channels were recorded: two total power 
channels (then averaged to a single total power Stokes 
I channel) and two correlation channels containing the 
Stokes Q and U parameters. We obtained 33 coverages. 

The telescope pointing was checked at time intervals 
of about 1.5 hours by making cross-scans of a nearby 
strong point source. The flux-density scale was calibrated 
by mapping the polarized source 3C138 and computing 
its total power flux densities of 2.49 Jy at 8.35 GHz using 
the formulae of Baars et al. I|1977[) . The data reduction 
was performed using the NOD2 data reduction package 
(Haslam I1974|l . The coverages were combined into final 
I, U, and Q maps using the spatial frequency weighting 
method (Emerson & Grave 1988). 

The distributions of Stokes parameters were combined 
into final maps of total and polarized intensity, polariza- 
tion degree, and polarization position angle. Noisy struc- 
tures smaller than the telescope beam-width were re- 
moved using a digital filtering. A convolution to the beam- 
width of 90" was applied to increase the sensitivity to 
the extended polarized emission. The r.m.s. noise level 
in the final map of polarized intensity at 8.35 GHz is 
0.039 mJy/b.a. 

NGC 4654 was observed, together with other large spi- 
rals of the Virgo Cluster, and a detailed analysis of their 
radio total and polarimetric properties will be discussed 
elsewhere. We also observed those galaxies at 4.85 GHz 
with a resolution of 2.5'. The contour map of polarized 
intensity of NGC 4654 at 8.35 GHz is presented in Fig.d 
with superimposed B-vectors of polarized intensity (i.e. 
E- vectors rotated by 90°), overlaid upon the Hi emission. 
The polarized emission shows strong North-South asym- 
metry that may impose a compression or stretching of 
magnetic fields on the southern side of the galaxy. Similar 
polarized structure is also visible in our lower-resolution 
4.85 GHz observations. 

The VLA observations represent part of a larger survey 
of 8 Virgo spiral galaxies that were observed in polarized 
radio continuum emission at 6 and 20 cm (PI: B. Vollmer). 
Since allof the data will be presented in a forthcoming ar- 

2 The 100-m telescope at Effelsberg is operated by the Max- 
Planck-Institut fur Radioastronomie (MPIfR) on behalf of the 
Max-Planck-Gesellschaft. 
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Fig. 1. Contours of polarized intensity of NGC 4654 at 
8.35 GHz with magnetic B-vectors proportional to the po- 
larized intensity superimposed onto the Hi emission. The 
contour levels are (3, 5, 8, 12) x 0.039 mJy/b.a. The map 
resolution is 90". Greyscale: Hi gas distribution (Phookun 
& MundyQlSJ). 

ticle, we only show the polarized radio continuum emission 
at 6 cm. 

NGC 4654 was observed at 4.85 GHz during 5 h on 
December 3 and 15, 2005 with the Very Large Array 
(VLA) of the National Radio Astronomy Observatory 
(NRAO) 3 D array configuration. The band pass was 
2 x 50 MHz. We used 3C286 as the flux calibrator and 
1254+116 as phase calibrator, which was observed every 
30 min. Maps of 17" resolution were made using the AIPS 
task IMAGR with ROBUST=3, which is close to pure nat- 
ural weighting (5). The maps were convolved to a beam 
size of the final maps of 20" x 20". We ended up with an 
rms of 10 /iJy/beam in total and polarized emission. In 
Fig. El the contour map of polarized intensity of NGC 4654 
at 4.85 GHz is presented, together with superimposed B- 
vectors of polarized intensity (i.e. E-vectors rotated by 
90°), overlaid upon the Hi emission. 

Similar to the Effelsberg data, the VLA data show a 
prominent ridge of polarized emission in the south of the 
galaxy. Moreover, this ridge extends to the southwest. The 
VLA data show that the lack of polarized emission in the 
Effelsberg data at the tip of the Hi head is due to strong 
beam depolarization since the magnetic field changes di- 
rection strongly following the Hi and optical envelopes 
(Fig. El- The magnetic vectors encircle the Hi head- like 
structure in the northwestern part of the galaxy. 

When we convolved the VLA data to the Effelsberg 
resolution, we found two major differences between the 

3 NRAO is a facility of National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 




Fig. 2. Contours of polarized intensity of NGC 4654 at 
4.85 GHz with magnetic B-vectors proportional to the po- 
larized intensity superimposed onto the Hi emission. The 
contour levels are (3, 5, 8, 12, 20) x 10 /itJy/b.a. The map 
resolution is 30". Greyscale: Hi gas distribution (Phookun 
& Mundv lTM5|) . 

VLA and Effelsberg maps: (i) the northern maximum 
is more extended in the convolved VLA map, which 
means that the Effelsberg telescope did not detect all the 
flux there (the northern maximum is smaller than the 
Effelsberg beam; cf. Fig.^), and (ii) the southwestern part 
of the main ridge of polarized emission in the Effelsberg 
data seems to extend into the extended Hi tail. Since this 
is not observed in the VLA data, this additional emission 
in the Effelsberg data is most probably due to the large 
beam. 

Both features: the asymmetry in the polarized emis- 
sion and the geometry of the magnetic field in NGC 4654 
should be explained by our MHD modeling and can consti- 
tute a powerful tool for distinguishing between numerical 
experiments and thus between different scenarios of gas 
dynamics (Vollmcr 2003) and the magnetic field evolution 
in this galaxy. 

3. Model 

The 3D MHD model of magnetic-field evolution (see 
Otmianowska-Mazur & Vollmer 120031 for more details) 
is applied to three selected numerical experiments of the 
Virgo galaxy NGC 4654 presented in Vollmer f2003) . We 
take only those cases into consideration where stars are 
treated as a non-collisional component (Vollmcr 2003): 
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Model GR: the gravitational (tidal) interaction with 
the second galaxy 
Model RPS: the ICM ram pressure. 
Model GRPS: the gravitational interaction and con- 
stant ICM ram pressure 

The Zeus3D MHD code (Stone & Norman ITMI i and b) 
is used to solve the induction equation 

3B 

— = V x (v x B) - V x (77 V x B) (1) 

where B is the magnetic induction, v the large-scale ve- 
locity of the gas, and 77 the coefficient of a magnetic diffu- 
sion (see Otmianowska-Mazur & Vollmer 2003). Realistic, 
time-dependent gas-velocity fields are provided by the 
3D N-body sticky-particle computations of Hi cloud com- 
plexes evolving in a galactic disk, a bulge and a halo sim- 
ulated by non-collisional particles (see Vollmer 2003 and 
Otmianowska-Mazur & Vollmer 2003). The clouds collide 
inelastically, and in the two experiments, they are affected 
by the ICM ram pressure due to motion of the galaxy with 
respect to the ICM. 

The velocity obtained from the N-body code has a dis- 
crete distribution. In order to solve the induction equation 
(Eq. ^) with the help of Zeus3D code, we have to interpo- 
late these velocities to the regular grid. In Otmianowska- 
Mazur & Vollmer (2003) we applied a spline function with 
a density-dependent smoothing length. It turned out that 
we had to use a large smoothing length to suppress the 
noise in the velocity field of the outer disk, which was due 
to a small, local particle density. In this way we avoided 
numerical artifacts of the magnetic field distribution in 
the outer disk layers, but at the same time we obtained a 
magnetic field distribution beyond the edges of the gas dis- 
tribution. For these reasons we decided to apply another 
method to solve the interpolation problem. The method 
is based on the kriging interpolation (see e.g Isaaks & 
Srivastava ll989(l modified to be used effectively in our 3D 
case. 

3.1. Kriging interpolation and our modification 

Kriging is a family of linear least-square estimation algo- 
rithms. It is used to approximate or interpolate data. The 
end result of kriging is to obtain the conditional expec- 
tation as a best estimate for all unsampled locations in 
a field and, consequently, a minimized error variance at 
each location. The kriging interpolation assumes (see e.g 
Isaaks & Srivastava 119891) that the value of any (vector) 
field at a point y(x) can be approximated by a weighted 
linear combination of a discrete set of vectors given at a 
certain set of points yi = y(xi) (provided by the N-body 
simulations, in our case) 

N 

V[x) = ^2viWi(x) (2) 

i=l 

where the weight coefficients Wi (x) determine to what ex- 
tent yi influences y(x). The best approximation condition 
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requires that the weights' sum has to be normalized by 

N 

Y] Wi(x) = 1. They can be determined from the set of 

i=l 

equations 

N 

*^2dijWi(x) + X(x) = dj(x) (3) 

where = d(Xi,Xj) and di(x) = d(xi,x) are called 
"covariance coefficients". These coefficients describe how 
closeby the points x% and Xj (or Xi and a;) are distributed. 
In our 3D case, we chose d falling exponentially with the 
distance dy = e~^ Xi ^ Xj ^ and di(x) = e~^ Xi ~ x ^. A slack 
variable (the Lagrange coefficient) X(x) is introduced here 
to balance numbers of unknowns with N + 1 equations 
(N for each Wi(x) plus the normalization equation). This 
method implies determining (to solve a set of N + 1 equa- 
tions) Wi(x) and X(x) for each point of interpolation x — 
in 171x171x71 points in our case. 

We noticed that we can substantially reduce the cal- 
culation time. If we define a new set of unknowns g^i = 
0, 1, . . . , N — 1 and a new slack variable k fulfilling the set 
of equations 

jV-l 

d ij9i + « = Vj (4) 

i=0 

JV-l 

E 9i = 0, (5) 

then, multiplying equations J3J) by gj and summing over 
j, we get 

jV-l N-l N-l N-l 

E Wil & E d ^ + A ( x ) E & = E to d A x ) ( 6 ) 

i=o j=o i=o j=o 

where we make use of the symmetry of dij. The second 
term vanishes (Eq. [SJ, and using Eq. (J3J), we can rewrite 
it as 

TV — 1 N-l 

E Wi(x) (yi - k) = ^2gjdj(x). (7) 

i=0 j=0 

Finally, by substituting Eq. J2J we obtain the following 
formula for the interpolated value 

N-l 

V( x ) = zl 9i d i( x ) + «• (8) 

Thus, for all interpolation points we solve set of equations 
(@J and (J5J) only once, and for each point x we calculate 
values di(x) and their scalar product with 

We applied this interpolation routine only to the 
velocity- field components. For interpolation of the mass 
distribution (used only for figures), we used a much 
faster method - smoothing particle masses with the 
Gaussian function with smoothing-length small enough to 
see gaseous features. 
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3.2. Model input parameters 

The dynamical models are fully described in Vollmer et al. 
(2003). As in our previous paper (Otmianowska-Mazur & 
Vollmer 2003), the induction equation is solved in rectan- 
gular coordinates (XYZ). We use 171x171x71 grid points 
along the X, Y and Z axis, respectively. The grid spac- 
ing is 400 pc in all directions, which gives the size of our 
modeled box as 68.4kpcx68.4kpcx28.4kpc. 

We assume the magnetic field to be partially coupled to 
the gas via the turbulent diffusion process (Otmianowska- 
Mazur et al. EODIi assuming the magnetic diffusion co- 
efficient to be 77 = 5 x 10 _25 cm 2 /s. The initial mag- 
netic field is purely toroidal with a strength of 10 /iG 
at a radius 2 kpc and falling linearly to zero toward 
R = 0. The magnetic-field strength is constant between 
2 kpc< R < 20 kpc and decreases with a Gaussian profile 
for 20 kpc< R < 30 kpc with a half width of 500 pc. It is 
zero for R > 30 kpc. In the vertical direction, the initial 
magnetic field strength is a Gaussian function with a half 
width of 500 pc. 

3.3. The polarization maps 

In order to obtain the simulated polarized intensity maps, 
we rotate the cube of 3D magnetic field according to the 
orientation of NGC 4654. Next, we integrate the transfer 
equations of synchrotron emissivities in Stokes I, Q, and 
U parameters along the line of sight: 



dl 



e 7 

Q I = I pei cos 2\ cos A — sin A 
U ) \ pei sin 2\ sin A cos A 




where the synchrotron emissivity is 



ej oc n re iB 



(7+l)/2 



(9) 



(10) 



We apply here the value of the relativistic electrons spec- 
tral index 7 = 2.8 and the Faraday rotation angle 



A oc A 2 / nth,B\\ds, 



(11) 



The subscripts _!_ and || are used with respect to the line- 
of-sight, and \ is the position angle of the sky-projected 
magnetic-field vector. The integration in Eq. ^] is along 
the line of sight, and observations wavelength is A. 

The intrinsic degree of synchrotron polarization is as- 
sumed to be p — 75%. Both thermal and relativistic elec- 
tron distributions are taken as Gaussian 



n re i\th °c e 



-(R/Ro) 



(12) 



where the radial scale-length Rq is set to 10 kpc and the 
scale-height zq to 1 kpc. We neglect Faraday effects by 
setting nth = 0. 

Finally for direct comparison with radio observations, 
the model maps of Stokes I, Q, and U are convolved with 
a 2-D Gaussian function with a HPBW of 20" and 80". 
The final Q and U maps are combined to obtain maps of 



polarized intensity and polarization angle (rotated by 90° 
to show the magnetic vector). 

As described in Sect. ^ the polarized radio-continuum 
emission traces the ordered large-scale magnetic field. In 
general, the total magnetic field of a spiral galaxy can 
be divided into two parts (see e.g. Beck et al. 1996): (i) 
a small-scale (smaller than a few hundred pc) magnetic 
field that is connected to the turbulent motions of the in- 
terstellar medium and (ii) a large-scale (larger than a few 
hundred pc) ordered regular magnetic field that is most 
probably due to a large-scale afi-dynamo. The large-scale 
magnetic field can be amplified by shear or compression 
motions of the interstellar medium where the magnetic 
field is frozen in. The small-scale magnetic field is ampli- 
fied by the turbulent motions of the interstellar medium 
(small-scale dynamo). 

The polarized radio-continuum emission is sensitive to 
a combination of magnetic-field compression or amplifi- 
cation. In addition, it is sensitive to the average field or- 
dering along the line of sight, which is a pure geometrical 
effect depending on the position and inclination angle of 
the galaxy. To separate these effects one can use the total 
power radio continuum map to calculate the degree of po- 
larization, which is only sensitive to the degree of intrinsic 
field ordering. 

As explained at the beginning of this section, we 
only model the large-scale magnetic field (Eq. . We do 
not have access to the small-scale magnetic field, which 
would require a model of the turbulent interstellar medium 
together with a small-scale dynamo. Therefore, we are 
not able to calculate a total power map for the model. 
However, we plan for the future to tie the distribution of 
the relativistic electrons (Eq. I12fl to the Ha emission dis- 
tribution. This emission traces ionizing O and B stars that 
will soon explode as supernovae giving rise to relativistic 
electrons. If we assume that the star formation history is 
approximately constant over the past ten million years, 
the Ha emission represents a reasonable measure for the 
relativistic electron distribution. We can try to link the 
total power map to the Ha emission using equipartition 
between the turbulent kinetic energy of the interstellar 
medium (which is connected to the star formation rate 
and thus to the Ha emission) and the total magnetic en- 
ergy. 

4. Results 

For the moment, we use the polarized radio-continuum 
emission for the comparison between our models and ob- 
servations. We show in the following that the information 
contained in the polarized radio-continuum is enough to 
distinguish between the different dynamical models. 

4.1. Gravitational interaction alone (model GR) 

Figure [3] shows the evolution of the polarized intensity 
distribution for the model of tidally interacting galaxies. 
The timcstcp of each snapshot is plotted in the upper 
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Fig. 3. Evolution of the polarized intensity at chosen time steps for the model only with gravitational interaction 
(model GR). Contours: polarized intensity in logarithmic scale. The polarized intensity vectors are superimposed onto 
the gas density distribution in logarithmic scale (grey plot). The contours are plotted at levels of (1, 3, 10, 30, 100, 
300) (in arbitrary units) in all plots. 

left corner of each panel. The closest encounter between steps. Figure 12 of Vollmer ( 2003} also shows the posi- 
NGC 4654 and NGC 4639 occurs at t = 460 Myr. The im- tions of the interacting galaxy NGC 4639. We use the 
ages are convolved to a resolution of 20" HPBW. This fig- inclination and position angle of NGC 4654 (i = 50° , 
ure can be directly compared to Fig. 12 of Vollmer (2003), PA — 120°). 
because we present snapshots exactly at the same time 
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The magnetic field distribution starts with an axisym- 
mctric, purely azimuthal configuration. Since the polar- 
ized radio-continuum emission is subject to Faraday ro- 
tation, i.e. the rotation of the magnetic field vector in 
the plane of the sky due to the magnetic field compo- 
nent along the line of sight, projection effects are im- 
portant. The magnetic field component along the line of 
sight is largest at the two ends of the galaxy's major axis. 
Therefore, Faraday depolarization suppresses the polar- 
ized radio-continuum emission there. Since the polarized 
radio-continuum emission depends on the line-of-sight in- 
tegration of the component of the magnetic field vector 
that lies in the plane of the sky, geometrical effects are 
important. And since this component of the magnetic- 
field vector is changing direction rapidly along the major 
axis, the corresponding polarized radio-continuum emis- 
sion is suppressed there. As a result, the synchrotron emis- 
sion distribution shows two maxima to the northeast and 
southwest of the major axis. The first snapshot in Fig. 
shows the polarized radio-continuum emission calculated 
using a slightly evolved magnetic field (t = 300 Myr). 
During the galaxy evolution, non-azimuthal gas flows in- 
duce radial components of the magnetic field resulting in 
the non-zero magnetic pitch angles without applying an 
explicit a dynamo mechanism. This effect becomes visible 
as a displacement of the polarized intensity maxima along 
the major axis (see e.g. Knapik et al. 2000). 

At t > 460 Myr, the gravitational interaction with 
NGC 4639 perturbs the gas flow leading to regions of 
gas compression and/or shear. Due to electromagnetic in- 
duction (Eq. ^) , the magnetic field is enhanced in these 
regions. Thus, the polarized intensity maxima follow gas 
compression and shear regions, but because the timescale 
of turbulent diffusion is larger than the viscous timescale 
of the gas, there is no one-to-one correspondence between 
them. In some places and at some timesteps, gaseous and 
polarization features coincide, while in other places and at 
other timesteps, arms of polarized radio continuum emis- 
sion lie just between gaseous ones (e.g. the northwestern 
maximum at t —700 Myr). More detailed analysis of po- 
larized intensity distribution shows that in general the 
maxima are left behind gas compression regions. This is 
not clearly visible in polarized intensity plots, as details 
are smoothed-out due to the line-of-sight integration and 
beam convolution. 



Further evolution (t > 800 Myr) shows a weakening of 
the magnetic field and thus a weakening of the polarized 
radio-continuum emission. In regions where gas falls back 
onto the galactic disk such as at the basis of the southern 
tidal tail at t > 900 My, the magnetic field is efficiently 
amplified, and maxima of polarized emission reappear. In 
the final time step (t = 1000 Myr), two comparable max- 
ima of polarized radio continuum emission coincide with 
gas compression or shear regions. 



4.2. Ram pressure alone (model RPS) 

The second simulation only includes a strong ram-pressure 
stripping event. The galaxy is diving deeply into the clus- 
ter potential where ram pressure becomes stronger due 
to the increased velocity of the galaxy and the higher in- 
tracluster medium density. The closest approach of the 
galaxy to the cluster center (M87), and thus maximum 
ram pressure occurs at t = Myr. The magnetic field 
evolution is strongly influenced by the gas compression 
and shear motions induced by ram pressure (Fig. HJ| . This 
figure can be directly compared to Fig. 3 presented in 
Vollmer (200S|), because we present model snapshots at 
the same time steps. The direction of the wind is iden- 
tical to Fig. 3 of Vollmer (2003). The inclination angle 
between the galactic and the orbital plane is i = 0°. The 
temporal ram pressure profile is a Lorentzian with a max- 
imum ram pressure of p ram = 5000 cm~ 3 (kms -1 ) 2 . In 
Otmianowska-Mazur & Vollmer ( 2003) we calculated the 
evolution of the distribution of polarized radio-continuum 
emission for a similar ram-pressure stripping event where 
the inclination angle between the galactic and the or- 
bital plane was i = 20° and maximum ram pressure was 
Pram = 2000 cm" 3 (km s~ 1 ) 2 . The evolution of the dis- 
tribution of polarized radio-continuum emission (Fig. |4*)l 
shows the same characteristics as that of Otmianowska- 
Mazur & Vollmer l)2003|l : due to the strength of the in- 
teraction, the magnetic field distribution and, in conse- 
quence, the maxima of polarized intensity follow almost 
perfectly gas compression regions for —50 Myr< t < 
200 My. Shortly before and after the maximum phase of 
the stripping event (t = —40 Myr and t = 80 Myr), we 
observe a bright maximum of polarized radio-continuum 
emission in the northeast of the galaxy center, at the lo- 
cation where the galaxy's ISM is compressed by ram pres- 
sure. At t = 200 Myr ram pressure and compression de- 
crease and rotation takes the magnetic field along with it. 
Since the galaxy rotates counter-clockwise, the result is an 
additional maximum of polarized emission in the south of 
the galaxy center. Rotation then makes the distribution of 
polarized emission approximately circular (t = 320 Myr). 
At t > 500 Myr the stripped gas that was not acceler- 
ated to the escape velocity begins to fall back onto the 
galactic disk in the southeast of the galaxy center (re- 
accretion). This leads to shear motion that enhances the 
large-scale magnetic field and thus the polarized radio- 
continuum emission. Re-accretion of stripped gas then be- 
comes stronger with time leading to a stronger maximum 
of polarized emission in the south of the galaxy center 
(see also Fig. 9 of Otmianowska-Mazur & Vollmer l2003f) . 
Differences between the polarized emission and the gas 
distributions are due to (i) the longer timescale of turbu- 
lent diffusion compared to the viscous timescale and (ii) 
orientation (integration along the line-of-sight) and beam 
depolarization effects. The final stage of the magnetic field 
evolution shows an approximately circular distribution of 
polarized radio-continuum emission with a maximum in 
the north and another maximum to the south of the galaxy 
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Fig. 4. Evolution of the polarized intensity at chosen time steps for the model with ram pressure event alone (model 
RPS). Contours: polarized intensity in logarithmic scale. The polarized intensity vectors are superimposed onto the 
gas density distribution in logarithmic scale (grey plot). The contours are plotted at levels of (1, 3, 10, 30, 100, 300) 
(in arbitrary units) in all plots. The arrow in the top left panel indicates the direction of the ram pressure wind, i.e. 
it is opposite to the galaxy's motion within the intracluster medium. 



center. The southern maximum is brighter than the north- 4.3. Gravitational interaction and ram pressure (model 
ern one. GRPS) 

Only when a weak, constant ram pressure is added to 
the gravitational interaction between NGC 4654 and 
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NGC 4639 are the observed Hi distribution and veloc- 
ity field of NGC 4654 reproduced fVollmer lM)3T> . When a 
small constant ram pressure is added, the magnetic field 
and thus the evolution of the polarized radio-continuum 
emission changes considerably (Fig.0. Again, it is possi- 
ble to compare Fig. directly to the snapshots of Fig. 14 
of Vollmer ((2003). The polarized intensity maps already 
show a clear asymmetry at early time steps because of gas 
compression due to ram pressure. Already at t = 300 Myr 
the polarized intensity maximum that corresponds to the 
southern maximum of Fig. [3] is weaker and displaced to 
the northwest. In addition, the polarized intensity maxi- 
mum that corresponds to the northern maximum of Fig. [3] 
is brighter due to gas compression and also displaced to 
the northwest. 

Around the timestep of closest approach between the 
two galaxies (t = 460 Myr), tidal forces exceed those due 
to ram pressure. Therefore, the distribution of the polar- 
ized radio-continuum emission for 450 Myr< t < 600 Myr 
is close to that of model GR (Fig. EJ). At later time steps, 
gas compression due to ram pressure in the northwest and 
gas re-accretion in the southeast dominate the gas dy- 
namics. This leads to a pronounced maximum of polarized 
emission in the northwest and an elongated ridge of polar- 
ized emission in the south of the galaxy center where the 
large low column-density Hi tail originates (see Phookun 
& MundyHES). 

4.4. Total magnetic energy evolution 

Figure shows the evolution of the total magnetic energy 
integrated over the whole cube normalized to its initial 
value for our three numerical models: GR (tidal interac- 
tion only), RPS (the ram pressure effect alone) and GRPS 
(tidal interaction and constant ram pressure). The mag- 
netic energy is calculated in a cylinder of 40 kpc diameter 
and 5 kpc height centered on the galaxy center. 

To investigate the impact of the initial conditions of 
the magnetic field on the evolution of the total magnetic 
energy, we used 3 different initial magnetic field configu- 
rations: 

(i) A "normal" disk: the magnetic field is constant up to 
a radius of 20 kpc, decreases with a Gaussian profile 
for 20 kpc< R < 30 kpc with a half width of 500 pc, 
and has a cutoff at 30 kpc. In the vertical direction 
the initial magnetic field strength is Gaussian function 
with a length scale of 500 pc and a cutoff at 1 kpc. 

(ii) A "thin" disk: the magnetic field is constant up to a 
radius of 15 kpc, decreases with a Gaussian profile for 
15 kpc< R < 20 kpc with a half width of 500 pc, 
and has a cutoff at 20 kpc. In the vertical direction 
the initial magnetic field strength is Gaussian function 
with a length scale of 300 pc and a cutoff at 600 pc. 

(iii) A "thick" disk: the magnetic field is constant up to 
a radius of 25 kpc, decreases with a Gaussian profile 
for 25 kpc< R < 40 kpc with a half width of 500 pc, 
and has a cutoff at 40 kpc. In the vertical direction the 



initial magnetic field strength is a Gaussian function 
with a length scale of 700 pc and a cutoff at 1.5 kpc. 

We observe that the evolution of the total magnetic en- 
ergy of the "normal" and "thin" disk initial conditions is 
very close. For t < 600 Myr the total magnetic energy 
of the "thick" disk model is systematically lower than 
that of the other models; i.e. the total magnetic energy 
grows more slowly than for the other two initial condi- 
tions. The reason is that the magnetic field amplification 
mainly takes place in the inner part of the galaxy's disk. 
Therefore, the ratio between the magnetic energy of the 
amplified field and that of the initial field is smaller for 
a larger "thick" disk than the same ratio for a smaller 
"thin" or "normal" disk. At t > 700 Myr, there is no 
significant difference between the total energies for the 
simulations with different initial conditions. We conclude 
that the evolution of the total magnetic energy depends 
only marginally on the initial configuration of the mag- 
netic field. 

The total magnetic energy initially grows in all three 
models (GR, RPS, GRPS). This effect can be explained 
by the magnetic field amplification due to gradients in 
the velocity field in early evolutionary stages. The initial 
field growth saturates at the time step of t ~ 450 Myr for 
models with a gravitational interaction (models GR and 
GRPS) and at t ~ 500 Myr for model with ram pressure 
alone (model RPS). These time steps correspond to the 
closest encounter with NGC 4639 (t = 460 Myr) and to 
the timestep of maximum ram pressure (t = 500 Myr). 
For t > 500 Myr the tidal interaction pulls gas out of the 
galaxy, taking the magnetic field with it beyond the lim- 
its of our computation box. Therefore, the total magnetic 
field energy decreases in models GR and GRPS. After the 
phase of strong interaction (t ~ 600 Myr), a part of the 
displaced gas falls back onto the galactic disk (phase of re- 
accretion, see Otmianowska-Mazur & Vollmer 2003) and 
its radial flows and shears connected to differential rota- 
tion of the galactic disk trigger magnetic field amplifica- 
tion very efficiently. As a result the magnetic energy of 
models GR and GRPS increases again at t = 600 Myr. 
At t > 700 Myr constant ram pressure modifies the gas 
dynamics creating shear and a compression region, which 
leads to a further enhancement of the total magnetic en- 
ergy, which then saturates at t = 750 Myr. 

Surprisingly, the enhancement of the total magnetic 
energy due to strong ram pressure compression (model 
RPS) is small (~ 10%). The characteristic timescale of 
the ram pressure stripping event (model RPS) is smaller 
than that of the gravitational interaction and re-accretion 
is stronger for model RPS. The subsequent shear mo- 
tion makes the large-scale magnetic field increase at t > 
600 Myr. With the diminishing of re-accretion and thus 
shear motions the magnetic energy decreases again at 
t > 800 Myr and stays approximately constant for t > 
1000 Myr. 

In all models we observe a significant growth in total 
magnetic energy at early interaction stages. Only for the 
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Fig. 5. Evolution of the polarized intensity at chosen time steps for the model with gravitational interaction and 
the constant ram pressure (model GRPS). Contours: polarized intensity in logarithmic scale. The polarized intensity 
vectors are superimposed onto the gas density distribution in logarithmic scale (grey plot) . The contours are plotted at 
levels of (1, 3, 10, 30, 100, 300) (in arbitrary units) in all plots. The arrow in the top left panel indicates the direction 
of the ram pressure wind; i.e. it is opposite to the galaxy's motion within the intracluster medium. 



model that only includes ram pressure (model RPS), the 
total magnetic energy at At = t — t > 300 (where to 
is the time of strongest interaction) is significantly higher 
(a factor of two) than the initial value (see Otmianowska- 



Mazur & Vollmer 2O0J5J) . The stripping of the large scale 
magnetic field, together with the ISM, might represent a 
way to amplify the intracluster magnetic field. The obser- 
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Fig. 6. Total magnetic energy evolution in our three mod- 
els - ram pressure alone (model RPS), only gravita- 
tional interaction (model GR), and gravitational interac- 
tion with constant ram pressure (model GRPS) 



5.1. The choice of the model polarized emission 
threshold 

As stated in Sect. 13.31 we can only compare the model 
to the observed polarized radio-continuum emission. The 
model polarized emission is calculated everywhere in the 
galaxy according to Eq. El For a proper comparison we 
need to know which level of the model polarized emis- 
sion corresponds to the observed r.m.s. noise level. Since 
we do not attempt to estimate the absolute synchrotron 
emissivity e/, because of our lack of knowledge about the 
absolute magnetic field strength and relativistic electron 
density, we decided to use a fixed fraction of the max- 
imum of the model polarized emission distribution as a 
model threshold (Fig. |SJ [SJ. To give an impression of the 
polarized radio-continuum distribution at lower flux den- 
sities, we show in Fig. the last snapshots of Fig. 01 0] 
andJSJwith one deeper contour level at 0.3 times the lowest 
contour of the corresponding snapshots. For the GR model 
(upper panel), the two maxima of Fig.|3]are connected by 
a ridge of faint, polarized emission. In addition, a faint 
emission region appears outside the main gas disk in the 
north, which is due to very low surface density gas (see 
Vollmer 2003). The faint component of the polarized emis- 
sion distribution of the RPS model (middle panel) forms a 
ring structure. Moreover, it extends a little bit farther out 
than the 3 times stronger emission. This is also the case 
for the faint emission of the GRPS model (lower panel). 

As a conclusion, the new polarized emission feature is 
due to very low surface-density gas in the case of the GR 
model. The lowest contours of the RPS and GRPS models 
(Fig. [TJ (i) do not show any new feature, (ii) are only a 
little more extended than the former lowest contours, and 
(iii) are a natural extension of the higher contour levels. 
Our method is thus not very sensitive to the choice of the 
model threshold. The same is true for the convolved data. 



vations of ICM show (see eg. Widrow 2002 and references 
therein) that such fields exist widely in clusters. 

We thus obtain a time-dependent global magnetic field 
amplification in all our models without any explicit mag- 
netic dynamo action (dynamo coefficient a = 0). The 
amplification was found to work very efficiently in the 
re-accretion phase - when the gas initially pulled out 
of the galaxy by an external interaction falls back due 
to the gravitational attraction by the host galaxy (see 
Otmianowska-Mazur & Vollmer 2003). 



5. Models versus observations 

We first compare our models to the low-resolution 
Effclsberg data. This comparison can already distinguish 
between the different models. The high-resolution VLA 
data confirm the choice of the "best fit" model. 



5.2. The Effelsberg data 

To compare our simulations to our Effelsberg observations 
(Fig. nj, we convolved the final model snapshots to the 
Effelsberg resolution of 80" HPBW at 8.35 GHz. These 
convolved snapshots can be seen in Fig. [S] together with 
our Effelsberg observations. All three numerical models 
lead to two maxima of polarized intensity, located north 
and south of the galaxy center. The second maximum of 
model RPS cannot be seen, because it is too weak. 

There are significant differences in the distribution 
and strength of the polarized intensity maxima among 
our models. The model only with gravitational interaction 
(model GR, Fig. |H1 top left) shows two maxima of almost 
equal intensity. The model with ram pressure alone (model 
RPS, Fig. |S1 top right) has only one bright maximum in 
the south and a very faint maximum in the north of the 
galaxy center that is 6.5 times weaker than the southern 
maximum. The simulation including a tidal interaction 
and a small constant ram pressure (model GRPS, Fig. [SI 
bottom left) shows two maxima with the northern one 1.3 
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Fig. 8. Models versus observations for the Effelsberg data. Contours: polarized intensity. Grey scale: gas distribution. 
The magnetic field vectors are shown as lines. Top left: model with a gravitational interaction alone (GR). Contour 
levels are (6, 8, 10, 12) x 2.5 (in arbitrary units). Top right: model with ram pressure stripping alone (RPS). Contour 
levels are (6, 8, 10, 12)xl5 (in arbitrary units). Bottom left: model with a tidal interaction and a small constant ram 
pressure (GRPS). Contour levels are (6, 8, 10, 12)xl0 (in arbitrary units). Bottom right: greyscale: Hi gas distribution 
(Phookun & Mundv lT995[l : contours: Effelsberg 8.35 GHz observations. Contour levels are (3, 5, 8, 12)x40 /iJy/beam. 



times weaker than the southern one. While the locations of 
the northern maximum arc the same for all three models, 
the location of the southern maximum is most eastward 
in model GR, westward in model GRPS, and in between 
in model RPS. 

5.3. The VLA data 

To compare the simulations to our VLA observations 
(Fig. 0), we convolved the final model snapshots to the 
VLA resolution of 20" HPBW at 4.85 GHz. These con- 
volved snapshots can be seen in Fig. Altogether with our 



VLA observations. Since the VLA image has a smaller 
r.m.s. noise level than the Effelsberg image, we plot deeper 
contours in the model polarized radio continuum emission 
maps. There is almost no difference between the low and 
high resolution maps of the GR model, i.e. there are two 
distinct maxima in the north and the south of the galaxy. 
The low-resolution map of the RPS model only shows the 
main ridge of polarized emission, whereas the high reso- 
lution map shows an additional faint ring-like structure. 
The difference between the high and low resolution model 
polarized emission maps is most significant for the GRPS 
model where the deep high-resolution map reveals a com- 
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Fig. 9. Models versus observations for the VLA data. Contours: polarized intensity. Grey scale: gas distribution. The 
magnetic field vectors are shown as lines. Top left: model with a gravitational interaction alone (GR). Contour levels 
are (3, 5, 8, 12, 20, 30)x4.3 (in arbitrary units). Top right: model with ram pressure stripping alone (RPS). Contour 
levels are (3, 5, 8, 12, 20, 30) x8 (in arbitrary units). Bottom left: model with a tidal interaction and a small constant 
ram pressure (GRPS). Contour levels are (3, 5, 8, 12, 20, 30) x8 (in arbitrary units). Bottom right: greyscale: Hi gas 
distribution (Phookun & Mundv H995[) : contours: VLA 4.85 GHz observations. Contour levels are (3, 5, 8, 12, 20, 
30) x 8 /xJy/beam. 



plex structure. The main ridge of polarized emission is 
located in the west and southwest of the galaxy center. In 
addition, a distinct maximum can be found in the north 
west at the edge of the gas disk. Polarized emission can 
be found within the entire southwestern quadrant of the 
galactic disk. The southern part of the polarized emission 
distribution splits into two arms in the southeast direc- 
tion. One arm is located within the inner gas disk, the 
other arm extends in the direction of the extended Hi tail. 



The inner arm is due to the tidal interaction, whereas the 
outer arm is mainly caused by ram pressure. 

6. Discussion 

The GR model produces neither the location of the south- 
ern maximum nor the intensity ratio between the two 
maxima. While the location of the southern maximum is 
reproduced best by model RPS, the intensity ratio be- 
tween the two maxima is much too large. Only the model 
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including a gravitational interaction with NGC 4639 and 
a constant ram pressure (model GRPS) reproduces our 
Effelsberg observations in a qualitatively satisfactory way. 
Both maxima are placed at about the correct positions, 
and the southern maximum is brighter than the northern 
one, although the intensity ratio is still too small (about 
1.3 in the model compared to 2.4 in our observations). 

Table 1. Mean magnetic pitch- angle values calculated in 
four chosen regions of the real galaxy and our three mod- 
els' disks (see Fig. |SJ) - ram pressure alone (Model RPS), 
gravitational interaction (Model GR), and gravitational 
interaction with constant ram pressure (model GRPS) 
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As a further step, we compare the magnetic field ori- 
entation in our models with that of our observations. We 
calculated the mean values of the pitch angles of the mag- 
netic vectors in four chosen regions of the galactic disk 
- one at the northern polarized intensity maximum, and 
three others in the region of the southern polarized inten- 
sity maximum: at the location of the maximum and on 
both sides (at a distance of one beam size, see Fig. 
These measurements are shown in Table ^ The models 
with a gravitational interaction alone (GR) and with ram 
pressure alone (RPS) yield pitch angles too small that are 
on the southern maximum. On the other hand, the pitch 
angle of the northern maximum of model GR is too large 
with respect to the observed value. The best agreement 
between the model and observed pitch angles for all posi- 
tions is found for model GRPS. 

The two prominent maxima of the high-resolution 
model's polarized emission map of the GR model are not 
present in the VLA observations. At first glance, the over- 
all polarized emission distribution of the RPS model seems 
to fit observations quite well. However, the faint northern 
maximum is not observed, there is no continuum emis- 
sion in the inner gas disk, and the southern maximum 
is about twice as strong as observed. Only the GRPS 
model's polarized emission distribution shows the key fea- 
tures of the VLA observations: (i) a distinct maximum 
can be found in the northwest at the edge of the gas disk. 
This maximum is due to gas compression by ram pres- 
sure; (ii) polarized emission can be found within the en- 
tire southwestern quadrant of the galactic disk; (iii) there 
is a region of strong polarized emission in the southwest of 
the galaxy. However, we do not observe the bifurcation of 
the polarized emission ridge to the southeast. Especially 
the arm pointing toward the extended Hi tail is missing 
in the observations. This might be due to a difference be- 
tween the model and real distribution of relativistic elec- 



trons. Despite this disagreement, we feel confident that the 
GRPS model reproduces our VLA observations best and 
thus confirms our findings based only on the Effelsberg 
data. 

We conclude that GRPS model not only successfully 
reproduces the Hi distribution and velocity field, but 
also successfully reproduces the observed polarized radio- 
continuum emission distribution. Observations of polar- 
ized emission can serve as an additional diagnostic tool for 
galactic interactions (see Otmianowska-Mazur & Vollmer 
2003 ) . It can add the information on gas flows perpendic- 
ular to the line of sight, adding a new dimension to radial 
velocity studies (see also Urbanik[2005). 

In the case of NGC 4654 we identified the northwest- 
ern polarized emission maximum as due to ram pressure. 
As in NGC 4522 (Vollmer et al. EES, we can detect the 
compression of the interstellar medium directly in polar- 
ized radio-continuum emission. This independently con- 
firms the claim of Phookun & Mundy (^95) that the 
linear Hi rotation curve in the northwestern part of the 
galaxy is due to compression and the same conclusion of 
Vollmer (2003) which was based on the comparison be- 
tween deep Hi data and a dynamical model. In the case 
of NGC 4522 (Vollmer et al. l!?Uo"6"jl . it was not even pos- 
sible to detect the compression region in the Hi velocity 
field. A part of the southern polarized radio continuum 
emission ridge is due to the tidal interaction, which is also 
responsible for the one-armed spiral. Another part is due 
to shear motions induced by the action of ram pressure. 
It is very difficult to determine regions of shear motions 
only on the basis of a radial velocity field. Since the po- 
larized radio continuum emission is sensitive to the 3D 
velocity field, one can gain important information about 
transversal motions. By definition, this information is not 
contained in a radial velocity field. Thus, polarized radio 
continuum emission gives us important information about 
shear and compression regions that may not be accessible 
via velocity fields. 

Our present study corroborates the results of Vollmer 
(2003), that NGC 4654 had a recent close and rapid 
tidal interaction with NGC 4639. Cepheid measurements 
(Gibson et al. l27)07)l Freedman et al. l2llrTT)) place NGC 4639 
several Mpc behind the Virgo cluster center (M87, see 
Sect. 1). However, H-band Tully-Fisher measurements 
(Gavazzi et al. H^9) yield a distance of NGC 4654 of 
1-2 Mpc in front of the Virgo cluster center. In light of 
our new results, we suggest that the Tully-Fisher distance 
is underestimated and that the line-of-sight distance of 
NGC 4654 is close to that of NGC 4639. 

Another result of our study is that a small amount of 
ram pressure is needed to reproduce the Hi gas distribu- 
tion, velocity field, and the distribution of the polarized 
radio-continuum emission. An intracluster medium (ICM) 
density of ~ 10~ 4 cm~ 3 is needed for this ram pressure 
(sec Sect. 1). If the ICM density distribution is spherical, 
the density at a distance of 5 Mpc is more than a magni- 
tude smaller. A density of 10~ 4 cm~ 3 is already reached at 
a distance of 1 Mpc from the cluster center (Schindler et 
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al. H999f) . We therefore suggest that the ICM distribution 
of the Virgo cluster is not spherical but elongated in the 
line of sight with an axis ratio of about 5:1 as observed for 
the galaxy distribution (Gavazzi et al. 1199911 Yasuda et al. 

7. Conclusions 

The Virgo spiral galaxy NGC 4654 was observed with the 
Effclsberg 100m telescope at 8.35 GHz and the VLA at 
4.85 GHz in polarization. We find an asymmetric distri- 
bution of the polarized intensity with a strong peak in the 
south of the galaxy. This shows that asymmetric polarized 
radio continuum distributions are typical of cluster spiral 
galaxies (Vollmer et al. 12001 Urbanik 12005) . 

The polarized radio continuum emission was used as 
an additional diagnostic tool for the interaction between 
NGC 4654 and its environment. To do so, we solved the 
induction equation for the dynamical model evolution of 
the interstellar medium of a cluster spiral galaxy and cal- 
culated the large-scale magnetic field in this way. With 
an assumed distribution of relativistic electrons, the dis- 
tribution of polarized radio-continuum emission can be 
calculated and directly compared to our observations. 
The input velocity fields are those from the simulations 
of Vollmer <|2003[) . who simulated three different inter- 
action scenarios: (i) a gravitational interaction between 
NGC 4654 and NGC 4639, (ii) ram pressure stripping by 
the intraclustcr medium, and (iii) a gravitation interac- 
tion, along with a small constant ram pressure. 

The direct comparison between our MHD simulations 
and the polarized radio-continuum observations shows 
that only the model with gravitation interaction and low, 
constant ram pressure can reproduce observations in a sat- 
isfactory way. We found that the low-resolution Effclsberg 
data are sufficient for distinguishing between the models. 
Thus, our study corroborates the results of Vollmer et al. 
( 200|J). This mixed interaction succeeds in reproducing the 
Hi gas distribution and velocity field (Vollmer 2003) and 
the polarized radio-continuum distribution. 

Since Cephcid measurements place NGC 4639 ~ 5 Mpc 
behind the center of the Virgo cluster, we suggest that the 
distribution of the Virgo intraclustcr medium is highly 
elongated thereby providing the necessary ram pressure. 
Thus, polarized radio-continuum emission, in combination 
with detailed MHD modeling, provides an important ad- 
ditional diagnostic tool for studying interactions between 
a galaxy and its cluster environment. This tool is comple- 
mentary to deep Hi observations. 
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RIGHT ASCENSION 

Fig. 7. Model polarized radio-continuum emission. Upper 
panel: GR, middle panel: RPS, lower panel: GRPS model. 
The contour levels are the same as in Fig. [31 El and 03 but 
the first contour is 0.3 times that of these figures. 



